Mt. Suhora M dwarf survey - Detection of eight short-period variable stars 



L. Fox Machado 3 *, A.S. Baran bc , M. Winiarski b , J. Krzesinski b , M. Drdzdz b 

" Observatorio Astronomico National, Institute) de Astronomia — Universidad National Autonoma de Mexico, Ap. P. 877, Ensenada, BC 22860, Mexico 
b Mt. Suhora Observaroty of the Pedagogical University, uL Podchorqiych 2, 30-084 Cracow, Poland 
c Iowa State University, 12 Physics Hall, Ames, IA 50011, USA 



Abstract 



o 

C 

(N 
(N 



The Mt. Suhora M dwarf survey searching for pulsations in low mass main sequence stars has acquired CCD photometry of 46 
M dwarf stars during the first year of the project (Baran et al 2011). As a by-product of this search hundreds field stars have been 
checked for variability. This paper presents our initial result of a search for periodic variables in field stars observed in the course 
of the survey. On the basis of the periodicity and the shape of the light curves, eight new variables has been detected, among which 
five are b Scuti stars and three likely RR Lyrae stars. Although variation in one of the stars has been previously detected, it was 
classified incorrectly. To support our classification, in August 2010, we performed spectroscopic observations to derive spectral 
types and luminosity classes for all eight variable stars. 
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TYC2594 1673 1, USNO-A2 120007693323, TYC 3527 2221 1, TYC 351766 1, USNO-A2 135002022405, BD+32 3331, 
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q 1. Introduction 

<J) , One of the advantages of using CCD cameras, compared to 
, photo-electric photometers, is multi-channeling, giving oppor- 
tunity to measure all stars in the CCD's field of view (hereafter 
i - 1 FOV). The number of stars may vary from a couple, in case 
of a small FOV and sparse fields, up to thousands for mosaic 
CCDs pointing near to the plane of the Milky Way. Appar- 
1^ ently, CCD technique gave us an opportunity to increase our ef- 
■^j- ficiency when an enormous number of stars has to be observed. 
This is particularly important in the sky surveys undertaken re- 
cently, like SDSS (York et al 2000), ASAS (Pojmariski 2001), 
and OGLE (Udalski et al. 1994). There are also surveys not 
aiming at covering all the sky, but focusing on specific fields or 
targets. A good example of these is the Kepler mission (Borucki 
et al. 2010). Some of the surveys are working in the time-series 
photometry mode, and these are the most suitable for variable 
star observations. 

Searching for new variables is also a secondary goal of the 
survey undertaken by Baran et al. (2011). Primarily, they are 
looking for stellar pulsations in M dwarfs with periods of about 
40 min. However, since the cadence period is very short (around 
10 sec on average), a flux variation with periods from a couple 
of minutes up to days, in adjacent field stars observed simul- 
taneously, can be easily detected. Here, we report a detection 
of seven new variable stars and one star with a wrong classifi- 
cation found in the literature. Amplitudes of these variability 
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are relatively large so the discovery were mainly done by eye- 
ball. Then, the Fourier technique (FT) and the Phase Dispersion 
Minimization (PDM) were applied to obtain periods of the de- 
tected variation. 

2. Photometric data 

Details on photometry taken in the course of the Mt. Suhora 
survey are given in Baran et al. (2011). To date, dur- 
ing the first two years of the survey, they observed sev- 
eral dozen M dwarfs with thousands of field stars. All field 
stars, up to reasonable brightnesses, which were different 
for each observing night, have been checked for variabil- 
ity. Eight stars, TYC 2644 8321, TYC 2594 1673 1, USNO- 
A2 120007693323, TYC 3527 2221 1, TYC 3517 661, USNO- 
A2 135002022405, BD+32 3331 and HD 188524 show bright- 
ness changes on a time scale of a few hours and were identified 
by eye, although USNO-A2 135002022405 and HD 188524 
show quite small amplitude variations. The list of the new 
variable stars detected during the survey is presented in Table 
1. TYC 35 17 66 1, was previously known to be a cepheid with 
a period of 1.585 d (Akerlof et al. 2000); however, its clas- 
sification was apparently wrong. Our higher signal to noise 
(S/N) photometry along with spectroscopy revealed that the 
light curve does not resemble that characteristic for Cepheids. 
Another star, HD 188524, was known to be a double star, 
but no photometry has been published in the literature. We 
have also searched the International Variable Star Index (VSXf] 
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database of the American Association of Variable Star Ob- 
servers (AAVSO) for known variables in our fields. None of 
them, but TYC 351766 1, has been assigned as a variable star. 

The light curves of all 8 variable stars are shown in Fig. 1 
through Fig. 6. In each panel we have indicated the name 
of the variable, the observing band, the date of observations 
(year-month-day, UT) and the observatory's name. Ampli- 
tudes of flux variations range from several parts per thousand 
(ppt = 0.921 mmag) for TYC26448321 and BD+323331, up 
to 40 ppt for TYC 35 17 66 1. Minimum to minimum periods, 
estimated by eye, range from less than an hour for USNO- 
A2 1 350 02022405, up to more than 5 h for TYC 35 17 66 1 . Pe- 
riods defined in this way are correct for pulsating stars but it is 
half of the period for binaries. Small amplitudes suggest that 
all stars reported here might be pulsators, although ellipsoidal 
variables or small amplitude binaries cannot be ruled out. Time- 
series spectroscopy could be very handy to distinguish between 
pulsators and binaries. 

Not all data have high S/N ratio. Stars presented here were 
not the target stars in the survey, so the exposure times have 
not been set to obtain high S/N for these stars. If their bright- 
nesses were comparable to those of the targets, we could expect 
good quality data. Otherwise, as in the cases of BD+32 3331, 
USNO-A2 135002022405 and HD 188524, a scatter is much 
higher than in other stars having comparable amplitudes of flux 
variations. More dedicated photometry of those three stars is 
needed for further analysis. 

3. Period analysis 

Since the data are sparse, the Fourier technique might not 
work very well, particularly for data with less than one period 
cycle. Therefore, we have also used Phase Dispersion Mini- 
mization (PDM) to better estimate periods. This tool is very 
useful for data sets with gaps or poor time coverage. In some 
cases PDM gives true as well as multiples of the true period. 
In those cases we picked the one close to our rough estimation 
which came directly from the light curve. We calculated peri- 
ods from 15 min up to slightly above the expected periods (by 
rough estimation) . 

The results of the PDM application is shown in Fig. 7. 
The vertical axe of the plots in this figure represents a vari- 
ance which is a indication of the goodness-of-the-fit. Any 
potential periodicities will occur as a minimum variance at a 
given period. We derived a 4.34 hour period for TYC 35 17 66 1 . 
The minimum variance value in this case is close to zero. 
It means that the data are of a very high S/N and there is 
one dominant period present in the data. For two other stars, 
TYC 2644 832 and TYC 2594 1673 1, the minimum variance 
we derive is also small, providing periods 2.80 and 4.35 h, re- 
spectively. TYC 3527 222 1 has two close minima of the vari- 
ance. It might be an indication of two close periodicities in the 
data. The deepest minimum gives 2. 1 1 h period. In the case 
of USNO-A2 120007693323 there are few local minima of the 
variance, with the deepest one representing the minimum-to- 
minimum variation. The period associated with this global min- 
imum is 4.19h. The Minimum value of the variance for three 



other stars with the poorest S/N ratios does not look that good 
as in the previous cases. This is caused by big scatter in the 
data, giving a variance value close to 1.00, even though the 
period estimation is still correct. That is a natural property 
of the PDM analysis. The periods for BD +32 3331, USNO- 
A2 135002022405 and HD 188524 are 0.88, 0.75 and 1.78h, 
respectively. 

4. Spectroscopic data 

We conducted spectroscopic observations at the 2.12 m tele- 
scope of the San Pedro Martir Observatory in Baja Califor- 
nia, Mexico during August 30, 2010 (UT). We employed the 
same equipment as in Fox Machado et al. (2010). In particu- 
lar, we used the Boiler & Chivens spectrograph installed in the 
Cassegrain focus of the telescope. The Thomson 2048 x 2048 
pixel CCD with a 0.14/im pixel size was attached to the spec- 
trograph. The 1200 lines/mm grating was used to cover a wave- 
length range from 3975 to 5300 A. We employed a dispersion 
of 0.6 A per pixel (44 A per mm) with a resolution of 2.6 A. For 
data reduction we followed the standard procedures using the 
IRAF package. 

Figure 8 shows normalized spectra of the seven variable stars 
discussed in the present paper. Because of technical difficulties 
we were not able to collect a spectrum for HD 188524. The 
derived spectral classifications are indicated in the caption of 
each panel. We marked characteristic spectral lines of these 
stars in each panel. As can be seen from Fig. 8, the spectra of 
TYC 2644 8321, TYC 35272221 1, TYC 2594 1673 1, USNO- 
A2 120007693323 and TYC 351766 1 show characteristics of 
F-type main sequence stars. A-type main sequence stars dif- 
fer from F-type main sequence stars because in the former the 
Balmer lines are remarkable for their strength, while in the later 
they are no longer pronounced. In addition, in F-type main se- 
quence stars the weak metallic lines are more numerous and 
stronger. 

A quantitative spectral assignment was made from the in- 
tensity ratio of calcium and hydrogen lines 4226 Ca 1:4341 Hy, 
since the 4226 Ca I line gets stronger toward later subtypes in 
F-type stars. We used the ratios of the following lines 4226 
Ca 1:448 lMg II and 4226 Ca 1:4385 Fe I, which show posi- 
tive luminosity effect, as a luminosity class indicator. Another 
commonly used feature in the spectral classification of the F- 
type stars is the G-band (lines i4305.6, i4308) - a molecu- 
lar band due to the diatomic molecule CH. It appears around 
type F3, depending upon the resolution of the spectrum, and 
strengthens toward later subtypes. As the figure shows, the 
G-band is present only in the spectrum of TYC 2594 1673 1, 
USNO-A2 120007693323 and TYC 35 17 66 1 (three panels on 
the right). 

USNO-A2 135002022505 shows a spectrum similar to that 
of TYC 2644 8321 -a F0V type star. However, the metallic 
lines are slightly weaker indicating an A9 type, instead. The 
spectrum of BD+323331, showed at the bottom panel of Fig. 8, 
is a typical A0V type star. It shows strong Balmer lines accom- 
panied with faint metallic ones. 
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An additional check of the spectral types was done by com- 
paring the normalized spectra with those of well classified stars 
available in the literature. In particular we used the library 
of stellar spectra STELIB by Le Borgne et al. (2003), whose 
spectral dispersion of 1 .7 A per pixel is about twice as low as 
ours. As an illustration, in the Figure 9 we compared the spectra 
of TYC 2644 832 1 (top) and STELIB FOV-type star HD 90277 
(bottom). The spectra are quite similar. 

5. Discussion 

Considering the spectral types, periods of the bright- 
ness variations and amplitudes of pulsation TYC 2644 8321, 
TYC35272221 1, BD+323331 HD 188524 and USNO- 
A2 135002022405 have been classified as 6 Scuti-type pul- 
sators . 

As well known, 6 Scuti-type variables are main-sequence or 
slightly evolved stars located in the lower part of the Cepheid 
instability strip in the Hertzprung-Russell diagram (H-R dia- 
gram). In general, the period range of 6 Scuti stars lies between 
0.5 to 8 hours and the spectral types range from A2 to F2. At 
present two main subclasses of 6 Scuti stars are recognized: 
the small-amplitude subclass and the high-amplitude subclass 
(Rodriguez & Breger 2000). The majority of the known 6 Scuti 
stars belong to the former subclass. 

The small amplitude 5 Scuti stars (SADS) are population I 
stars which pulsate with a number of nonradial p modes ex- 
cited to typical photometric amplitude of 0.02 mag, whereas 
the high amplitude 6 Scuti stars (HADS) are mainly Population 
II monotone radial pulsators with amplitudes bigger than ~ 0.3 
mag. While the HADS are low rotators, the SADS tend to have 
much larger vsini values, in some cases reaching 200 kms~' 
(Rodriguez & Breger 2000). Since the 6 Scuti-type stars are 
relatively bright objects, they can be easily detected from the 
ground. In fact, several 6 Scuti stars have been discovered ac- 
cidentally when taken as reference stars of observations of well 
known 6 Scuti stars [e.g. Fox Machado et al. 2002; Li et al. 
2002, Fox Machado et al. 2007]. 

In addition to the HADS and SADS there is another type of 
pulsating stars which fall in the same region of the H-R dia- 
gram as the 5 Scuti stars, namely the SX Phoenicis variable 
stars (SX Phe). The SX Phe stars are sub-dwarfs with a spec- 
tral type A2-F5. Their light variations resemble those of the 6 
Scuti variables, but with shorter periods from 0.5 h to 2 h and 
larger amplitude of pulsation, up to 0.7 mag in V. The vast ma- 
jority of all SX Phe variables are metal-poor stars and therefore 
have been found in globular clusters. The field SX Phe stars are 
rare, even thought are interesting targets for the study of radial 
pulsations (e.g. Fauvaud et al. 2010). 

TYC 2594 16731, TYC 3517661 and USNO- 
A2 120007693323 are ambiguous cases for classification. 
Regarding the light curve shape these variables seem to be 
monoperiodic radial pulsators, such as Cepheids, RR Lyrae 
stars, HADS and SX Phe variables. In contrast, their ampli- 
tudes suggest that all of them might be small amplitude 5 Scuti 
stars. Since their period is rather too short for a Cepheid star 
and too long for a SX Phe variable, these variability types can 



be ruled out. Moreover, they cannot be 5 Scuti stars by virtue 
of their late F spectral types. Therefore, they might either be 
RR Lyrae type stars or small amplitude binaries. In the former 
case, they could be RR Lyrae type stars of subclass RRc which, 
as well known, oscillate in the radial first overtone with a 
period in a range from about 0.2 to 0.5 days. In the later case, 
they could be W Ursae Majoris-type stars (W UMa) which 
are close contact binary systems. As it is known, the stars in 
W UMa systems are main-sequence stars of nearly the same 
spectral type, mostly from around middle F to early G. The 
W UMa systems typically have periods between 0.25 and 1.2 
days. 

We tentatively assigned them to the group of RRc-type vari- 
ables, but more dedicated observations are needed for a more 
convincing classification. 

6. Conclusion 

We have presented our results of a search for new variable 
stars in field stars observed in the course of the Mt. Suhora 
M dwarf survey project (Baran et al. 201 1). We have identified 
a total of eight variable stars. Seven of them are previously un- 
catalogued variables. To shed more light on the nature of these 
variable stars, spectroscopic observations were performed. 

TYC26448321, TYC 35272221 1, USNO- 

A2 135002022405, BD+323331 and HD 188524, have 
been classified as 5 Scuti-type variables taking into account 
not only the spectral types but also the oscillation periods 
and amplitudes of pulsation. For these stars photometric 
observations in a larger time span will probable reveal their 
multiperiodic behavior. 

USNO-A2 120007693323, TYC 2594 1673 1 and 
TYC 3517661 may be classified as either RRc variables 
or W UMa binaries. Spectroscopic time-series could help in 
disentangling their nature. 
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Table 1 : A list of new variable stars detected in the course of the Mt. Suhora survey. Only the star with the asterisk was known earlier but classified wrongly. 
Coordinates are given for the 2000.0 epoch. 



ID 


RA 
(J2000) 


DEC 
(J2000) 


V 

(mag) 


Sp. T. 


Period 
(hour) 


Class. 


TYC 2644 8321 


19:07:47 


32:35:02 


10.0 


F0V 


2.80 


6 Set 


TYC 2594 16731 


17:04:27 


32:19:46 


13.0 


F6IV 


4.35 


RRc 


USNO-A2 120007693323 


15:55:18 


35:18:19 


12.0 


F6IV 


4.19 


RRc 


TYC 3527 2221 1 


18:35:40 


45:40:15 


9.9 


F0V 


2.11 


6 Set 


*TYC3517661 


17:36:28 


49:05:20 


11.0 


F9V 


4.34 


RRc 


USNO-A2 135002022405 


02:09:14 


49:27:18 


13.0 


A9V 


0.75 


6 Set 


BD+32 3331 


19:08:24 


32:29 13 


9.6 


A0V 


0.88 


6 Set 


HD 188524 


19:53:42 


44:36 


10.4 


F0V 


1.78 


6 Set 



TYC 2644 832 1, R, 2010-06-04, SPM 




1352.80 1352.85 1352.90 1352.95 1353.00 




1353.80 1353.85 1353.90 1353.95 1354.00 



TYC 2644 832 1, R, 2010-07-01, Suhora 




1379.40 1379.45 1379.50 1379.55 



TYC 2644 832 1, R, 2010-07-03, Suhora 




1381.40 1381.45 1381.50 1381.55 



TYC 2644 832 1, R, 2010-07-12, Suhora 

u 30 E ^ ; — ^ ; 




1390.40 1390.45 1390.50 1390.55 

HJD - 2454000.0 



Figure 1 : Light curves of TYC 2644 832 1 . Data have been taken on five nights. Note different S/N ratios for first two and for the rest of the nights as a reflection of 
the two different telescopes used. The captions in this and the following figures contain the name of the star, the filter used, evening date of observations, and site: 
SPM - San Pedro Martir Observatory, Suhora - Mt. Suhora Observatory 
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TYC 2594 '673 1, R, 201 0-06-28, Fick 



-20^ 

1376.60 



s2w c 00 c 



1376.65 



1376.70 



Sf§>% 



•'in; 



Cq C 



1376.75 



1376.80 



1376.85 



1376.90 



20 



TYC 2594 '673 1, R, 201 0-07-02, Fick 
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Figure 2: Same as in Fig.l but for TYC 2594 1673 1. Fick Observatory. Data have been taken on two nights. 
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USN0-A2 1200 07693323, R, 2010-07-01, Fick 
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USN0-A2 1200 07693323, R. 2010-07-02, Fick 
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Figure 3: Same as in Fig.2 but for USNO-A2 120007693323 
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Figure 4: Light curves of three variables detected TYC 3527 2221 1, TYC 3517 66 1, and USNO-A2 135002022405. TYC 3517 661 was previously classified as 
cepheid. 
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BD+32 3331, R, 2010-07-01, Suhora 
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Figure 5: Light curves of BD +32 3331. Data have been taken on three nights. 
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HD 188524, R, 2010-08-11, Suhora 
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Figure 6: Light curves of HD 188524.This star has the worst S/N ratio among all stars presented in this paper because the target star was significantly brighter. 
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Period [sec] Period [sec] 

Figure 7: PDM for all eight variable stars detected. In some cases we included only one night, even though we collected data on more nights. 
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Figure 8: Normalized spectra of the seven variable stars. We provide a spectral classification in each panel. 
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Figure 9: A comparison of the TYC 2644 832 1 spectrum with that of a FOV-type star from the STELIB library (see text for details) in the wavelength region between 
3950 -5250 A. 
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